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Abstract 

The system consisting of a self gravitating perfect fluid and scalar field is consid- 
ered in detail. The scalar fields considered are the quintessence and "tachyonic" 
forms which have important application in cosmology. Mathematical properties 
of the general system of equations are studied including the algebraic and dif- 
ferential identities as well as the eigenvalue structure. The Cauchy problem for 
both quintessence and the tachyon is presented. We discuss the initial constraint 
equations which must be satisfied by the initial data. A Cauchy evolution scheme 
is presented in the form of a Taylor series about the Cauchy surface. Finally, a 
simple numerical example is provided to illustrate this scheme. 
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1 Introduction 



In the arena of cosmology and stellar structure, the importance of the perfect fluid 
source cannot be understated. The assumption of large-scale isotropy will demand, 
via the field equations, that the stress-energy tensor supporting the universe possess 
the algebraic structure of a perfect fluid. As well, cosmologies containing scalar fields 
have long been considered. The motivation being that the scalar field would represent 
some exotic component of matter, which could explain various puzzling phenomena 
presented by observational cosmology. More recently, viable alternatives to the stan- 
dard big bang theory have been put forward, the most promising of which involve an 
inflationary era at some early time. Most of these theories invoke a scalar field to 
play the role of the inflaton in the early universe. At late times, the scalar field may 
again play an important role as a "dark energy" field. It has been observed that the 
universe has recently entered an acceleration phase and some exotic dark energy must 
presently dominate PP, It is well known that a scalar field provides a simple model 
to explain this acceleration. 

The scalar field may take many forms. Aside from the traditional quintessence 
field, which arises from a Lagrangian motivated by relativistic continuum mechanics, 
there is also the recently popular "tachyonic" scalar field from string theory, which 
possesses a Born-Infeld type action. Both of these fields have been used extensively 
in cosmology as the inflaton, dark matter and dark-energy (see |3], [H], [7j, 
0, 0, Cni, CH, Cni, mi, ini, Cni, im in and references therein). Usual 
cosmological models are spherically symmetric. However, more general cosmological 
models, admitting a three-parameter group of motions (Bianchi type I-IX) have also 
been studied. 

For the above reasons, we believe that it is important to study the combined 
system of Einstein-scalar field-perfect fluid equations in detail. Specifically, in section 
2 we present the system in its full generality for both the quintessence field and tachyon 
field. We discuss the system of equations and the number of functions which can be 
apriori prescribed. This is dictated by the number of functions present as well as 
the number of equations and identities. The eigenvector and eigenvalue structure is 
briefly discussed as well as its relation to physical quantities such as energy density 
and pressure. 

The Cauchy problem has long been studied in general relativity and is of great 
use due to the difficulty involved in obtaining exact solutions to physical problems 
within the theory (see jTHI, [20]; as well as where the gravitating complex scalar- 
Maxwell field system has been studied). In solving Einstein's equations, one often 
considers physically acceptable initial data (for example: data from plasma physics 
for stellar collapse, a specific form for the scalar potential inspired by string theory or 
particle physics in cosmology, etc.) The future evolution of this prescribed data is to 
be determined, ideally by analytically solving the Einstein equations and matter field 
equations of motion. However, this is often not possible and either some simplifying 
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assumptions must be made or else approximation methods must be utilized. It is 
useful therefore to numerically evolve the initial data of the full, unsimplified model 
to determine the behaviour of the system at later times and, in some cases, to aid 
in determining what form the exact solution should take. A thorough treatment 
of the system consisting of a massless self-gravitating scalar field in the context of 
gravitational collapse (in the retarded time gauge) may be found in [22]. 

The Cauchy problem in general relativity is non-trivial. Initial data is not 
freely prescribable due to the existence of four constraining equations for the ini- 
tial data. (These arise due to the differential identities among the system of field 
equations.) The constraint equations for both the quintessence plus perfect fluid and 
tachyon plus perfect fluid systems are derived in this paper along with the functions 
which one, in principle, may prescribe. In section 3 the Cauchy problem for the 
Einstein-perfect fluid-scalar field is properly posed using geodesic normal coordinates 
(which are particularly well suited to cosmological studies). A scheme is presented 
which will evolve the this data from the initial t = tQ = constant hyper-surface ana- 
lytically into the bulk {t > to)- Finally, we end with a simple example to illustrate 
the use of the scheme and its convergence properties. For an excellent review of the 
Cauchy problem in general relativity, the reader is referred to |23j . 



2 General theory 
2.1 Quintessence field 

The system considered consists of gravity coupled to a perfect fluid and quintessence 
scalar field, (j). The fluid energy density, pressure and scalar potential are denoted by 
p, p and V{(j)) respectively. Derivatives with respect to the scalar field are denoted by 
a prime. 

The quintessence Lagrangian density has the form ^: 



(1) 



^Conventions here follow that of ^3] 
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The Einstein and supplemental equations for this system are given by: 



S^, := G^, - MTt = 0, 



K + {P + P) u^u^^v^%, 



T ■=T^ =0 
■^y ■ u]ij, 

a := <p-f^ - V'{<P) = 0, 

U := u'^u^ + 1 = 0, 

/C:=[(p + p)n"]^„-ti>,„ = 0, 

■= (P + P)u°'u^;a + {SI + U°'Uu) P,a = 0, 
C{ga/3, dap a) = 0- 



(2a) 

(2b) 

(2c) 
(2d) 
(2e) 
(2f) 

(2g) 
(2h) 



(The denote four possible coordinate conditions.) The algebraic and differential 
identities are: 



+ SvrT^ = 0, 
% - {a(f)^^ + u„IC + Tu) = 
1 
2 



(p + p) {U - 1)^^ + Uu'p,, - u'T. = 0. 



In this system of equations the number of unknown functions is 

10 (g^p) + 1 (0) + 1 {V{<P)) + 1 (p) + 1 (p) + 4 {un = U 



(3a) 
(3b) 

(3c) 
(4) 



whereas, from ()2all2hp . there exist twenty-five equations. However, there are nine 
identities fl3all3cp and therefore only sixteen independent equations yielding an under- 
determined system. We can therefore prescribe two functions out of the eighteen. In 
case an equation of state is imposed. 



S{p,p) 



dSip,p) 
dp 



+ 



dSjp, p) 
dp 



= 0, 
> 0, 



(5) 



we can still prescribe one function. 

We can now explore eigenvalues of T^. Assuming that 

0,^ ^ 0, 
u°s^ = 0. 



(6) 
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We obtain from (|2ljl Oe|) and (0) that 



u 



1 



(7a) 
(7b) 
(7c) 



It is evident from the above equations that, in general, none of the vectors -u^, 0'^, 
are eigenvectors. However, consider two special cases. In case-I, we take the scalar 
gradient perpendicular to the fluid velocity and s^: 

u''(p,u = 0, (8) 
s'd. = 0. 



By (fraTEl and 



we get: 



7^" o/3 



(9a) 
(9b) 
(9c) 



In this case, m^, 0^ are eigenvectors and is a two-fold degenerate eigenvector. 
The corresponding eigenvalues in ()9all9c|) represent the proper mass density and two 
principal pressures, respectively, of the matter. A positive value of the potential, 1/(0), 
make a positive contribution to the energy density and negative contribution to the 
pressures. The structure here is algebraically similar to the anisotropic fluid. 

In case-II, assume that the scalar gradient is colinear with the fluid velocity: 

0'^ = bu^", (10) 
b ^ 0, = 0. 

By equations ()7all7c|) and (fTn|) . we arrive at: 



r^0^^ 



1 



P - \<l>'^<t^,p + v{4>) 



p--0^^0,;3-V(0) 



U 



(11a) 
(lib) 
(11c) 



In this case, m^, 0^ and are all eigenvectors. The structure here is algebraically 



similar to a perfect fluid with proper mass density p 



i0^^0,;3-y(0). 



+ y(0) and pressure 



2.2 Tachyonic scalar field 

The tachyon field Lagrangian density is given by: 



with $ representing the tachyon field. 

In the case of the tachyonic scalar field, the governing equations read: 



:= G"t - = 0, 



T ■=T^ = n 



^1 + 



+ {p + p) u^'uu+p6''^, 



a :- 



y'($)^/rT¥^ = o, 



:= + 1=0, 



(12) 

(13a) 
(13b) 
(13c) 
(13d) 

(13e) 
(13f) 

(13g) 
(13h) 



Again the system is underdetermined by two which allows us to prescribe an equation 
of state and one quantity (usually p or p). Defining vectors as in ^ we find: 



=y($)- 



1 f^S 



p- 



1 + 



i^($)v/i + + P 



y($) 



+ L _ F($)^l + $;«$,^ 



(14a) 
(14b) 
(14c) 



As with the quintessence field, none of the vectors tt^, </)'^, are eigenvectors. In case 
the vectors are mutually orthogonal the following equations hold: 



T%uf = - [p + y ($) yrT¥^ 



P 



1 + 



p-F($)v/i + w,. 



(15a) 
(15b) 
(15c) 



Note that in this case, not only are u^, (j)'^ and eigenvectors, but the stress-energy 
tensor possesses the algebraic structure of an anisotropic fluid. 
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In case the fluid velocity and scalar gradient are colinear, the stress-energy 
tensor has following structure: 



p + 
p + 



^1 + $;^$_^ 



1 pb 



which is similar to that of a perfect fluid. 



(16a) 

(16b) 
(16c) 



3 The Cauchy problem 

Here we present the Cauchy problem for the self-gravitating system of a perfect fluid 
and scalar field. 

Let us consider a contractible space-time domain (see figure 



D 

X 



Dx (ti,to) 
(x,t) G D. 



(17) 



(D, t.) 




Figure 1: Domain for the Cauchy problem. Initial data is specified on the hyper-surface 
domain {D,to) and is evolved to the hyper-surface domain {D.,ti)- 



Let a differentiable symmetric tensor field S^v{x) exist in D of the space-time. We 
cite Synge's lemma 
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Lemma 3.1 Synge's lemma: Let S^^^{x) he a symmetric, differentiable tensor field 
in the domain D of space-time. Then, the following two statements are mathematically 
equivalent: 



A: S^i,{x) = for X e D . 
B: Sij{x) - ]^gij{x)S1^{x) = 

and S"^.^ = in D, with 5''^(x, to) = 0. 



(18) 



The above lemma can be applied to Einstein's field equations. In that case, the 
following two statements are mathematically equaivalent: 



A: 
B: 



in D 
1 



^ijix) - -gij{x)Sl{x) = 
and 7^ = in D, with S^^{-x,tQ) 



(19) 



0. 



We shall investigate equations (fT^ involving a quintessence scalar field and perfect 
fluid in the following geodesic normal coordinates: 



ds^ - 



J- jk 



= 9tiu{x) dx^dx^ = — ((ix°)^ + gij{x) dx^dx\ 
--■ 9ij{^), 



^ jk^ ^ Ok 



I to 



)9 9ik,oi ^ i'l 



5 9ij,o,o\^) '■- 



1 



9ij,o,o{^j t) 



I to 



(20) 



(21) 



To clarify, glj will subsequently be used to denote the bulk spatial metric and g^j will 
be used to denote the spatial metric on the initial hyper- surface. 

Computing the Ricci and Einstein tensor components from ()19p we obtain: 

R,,{x) =Rl^{^, t) + \gl,^, + -J''gl,gl, - \9^''9lA^. (22a) 



Ginix] 



--Q 



m 



' i _ t 

9kl,0:ii 9iin.0 



9 9kL0 



2 1 

_ _Jkm Jin J J 

qU y 9kL09mn,0- 



(22b) 

(22c) 



Here, sharp symbols (jj) on indices denotes covariant differentiation with respect to 
the three dimensional metric g\y The Ricci and Einstein tensor components can be 
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expressed as (from ^ and dUi-IHcj)): 



Rij{^, to) =Rij{^) + ^^ij- o,o(x) + T/'^fc«,o^ii,o - l9''^9M,o9ij,o, (23a) 



1 



Goo(x, to) =2^(^) + 8 [t^9ki,o] - ^t""9^"9ki,o9mnfi- 



(23b) 
(23c) 



Note that the extrinsic curvature components of the initial hyper-surface are given by 



\gijQ and the above equations are also derivable from the Gauss-Codazzi 



equations [21, [23, ES- 



3.1 Quintessence field 

We shall now explore the field equations and subsequent evolution with quintessence 
scalar field and perfect fluid sources. Using Synge's lemma and consequent equations 
(Uni) and (j22il-l22cl) we arrive at: 

^ij ~ 7^9ij^ a = ^ij (^1 + 9 fi'ij.O.O + TS'" 9lmfi9ijfl ~ 7^9^ 9Hfi9mj,0 



- Stt + {p + p) UiUj + gl- 

1 



29^''^h<P,m + V{<P) + -{p-p) 



'g9 L 



- 0,0,0- (in v^)^ 0,0 -^'(0)=O, 



It 



?7(m) := v^TTT^X^ > 1, 



0, (24a) 
(24b) 
(24c) 



/C = f/(«)p,o + ixV,^ + { (V?^') ^ + (v^) ^ + V^f/,o} = 0, (24d) 

^0 = [1 - iU{u)f] p,o - u'U{u)p^k -ip + p) !^u% + ^g^'"^gL,ou''u,~^ - UUo = 0, 

(24e) 

= (p + p)f/(M)M,, - ^(P + pWg^^'^gl^^oUb + [p,, + + Uu,p,o\ = 0. (24f) 

Moreover, on the initial hyper-surface {t = to), the initial data must satisfy: 
1 



£io(x,to) = -^^""(x) [g^^^o-ri - 9ni,0;rn\ ' [0,i0,O - {p + p)UUi\^^^ , 
foo(x,to) =Ji?(x) + i [r"(x)^^„,o]' - Ir'^r^gmnfi-gabfl 



(25a) 



8vr|^[(0,o)' + /7"''V,a0,6]+m + (P + P)t^'-p}^ =0. (25b) 
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The initial data is to be prescribed as: 

p(x,to) =/u(x), p{x,to) = r]{x.), (26a) 

Ma(x,to) =«'a(x), (26b) 

0(x,to) =x(x), 0,o(x,t)|to = ^(x), (26c) 

9abi^) =9ibi^, to) = 7afe(x), (26d) 

9ib,oi^,t)\to='^ab{^)- (26e) 

The given functions /i, a, Wa and ^ are of class in D_ G M.^ , the functions "f/^ab are of 
class and the functions jab are of class C^. Moreover, the consistency conditions 
al- l25b|) must be satisfied as: 

^7™" [i^mn;a ' ^na;Tn] + StT {^X,a ' [^^ + V] U Wa] = 0, (27a) 



— »7r 



^ (e' + r\,ax,b) + v^(x) + (/i + ?7) - 77 



0. (27b) 



(Here, barred quantities are derived from the prescribed metric 7mn(x) ). The system 
of three dimensional partial differential equations is underdetermined so that infinitely 
many solutions should exist locally. Restricting the field equations to the initial hyper- 
surface, t = we can rearrange the equations symbolically as: 

-^jj Q = initial data and their spatial derivatives, 

</',0,0 1*0 — ) 
f^P,0|to = ' 

{{\-\j'\vA\,. = " " " " 

[fl + V]UUa,0\to = " " " " 

Thus, the higher derivatives of gij, (j), p, p, and Ua at the initial hyper-surface are 
determined by allowable Cauchy data and their spatial derivatives. Assuming the 
functions giji^.t), 0(x, t), p(x, t), p(x, t), and Ma(x, t) are real-analytic functions in 
to < t < t, we can determine arbitrary higher order derivatives with respect to 
by differentiating the field equations in ()24al - I24fjl and restricting subsequently those 
equations to the initial hyper- surface. Thus, the power series: 



fi'ii(x,t) 


=7ii(x) + ^^i(x)(^- 


to) + ^9ij,o,o{^-to? + - , 


(28a) 


0(x,t) 


=x(x)+e(x)(^-^o) 


+ ^(P, 0,0 \to{t -tof + ■■■ , 


(28b) 


p(x,t) 


=/i(x) + p,0\toit-to 


) + ••• , 


(28c) 


P(x,t) 


=V{^) +P,o\to{t - to 


) + ••• , 


(28d) 


Ua{x,t) 


=Wa(x) + Ua,0 \hXt - 


to) + ... , 


(28e) 
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can be generated. By the real analyticity conditions, there exist a ti > to such that 
all the power series in ()28all28e|) converge absolutely for all t G (ti,to] and converge 
uniformly for all t G [ti — 6, to] for a sufficiently small 6 > 0. 

3.2 Tachyonic field 

In the case of the tachyonic scalar field, the field and supplemental equations may be 
expressed as: 



- Stt <^ V{(t)) 



+ 4 (1 + 



+ (p + p)UiUj + -g\.{p - p) J> = 0, 



(29a) 



'9' 9 



1-1 



U 



- [1 + 
-2$,o<^>,i^""^<^>,n;o] = 0, 

=: U := ^/TTg^^h^b > 1, 



^-$,0,0- (lnv¥)^$,o-[ln|V^(0)|]' 

(^<l>,„^,, + ^i^o-^-o) ($,o)' $,0,0 



:= f/ p,o + MV,a + 



/ ,a 



0, 



(29b) 
(29c) 

(29d) 
(29e) 



1 



= (p + p)t/(ti)ti,,o - ^(P + P)t^^7"""i,,o«f> + [P,j + u\P,b + ^«,P,o] = 0. (29f) 



As before, there exist constraints on the initial data on the t = to hyper-surface: 

^ao(x, to) = "~ -^9 [9mn,0]a ~ 9na,0;m\ 

<^',a$,0 



- Svr <^ y(0) 



1 + c/tt'^^$,,$,, - ($,, 



11 2 1 

^Oo(x, to) ^2'^''' 8 [^"""^"^"'O] ~ g9'™'"fi'" 9mn,09ab,0 



^TT <^ y(0) 



1 + ^«"^<l>,a<l>,fc 



-ip + p)Uua} =0, (30a) 

|io 



+ (p + p)?72 -p j> =0. 

(30b) 
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The initial data is prescribed as in the quintessence case ()26a|) - ()26e|) with replaced 
by $ (differentiability requirements remain the same). In this case, the consistency 
equations (|3(Ja|) - ()30b|) become: 



+ r])U Wa} =0, 

(31a) 



ab 



+ S7f { V{x) 



1 + 7°''X,aX,fe 



y 1 + 7™"x,™x,n - e 

At this point, the evolution is governed by the equations (j28aj) - (j28ep as before. 



(31b) 



3.3 A specific example 

It is instructive to demonstrate how the above scheme works with an explicit example. 
We shall study the numerical evolution of a system whose analytic properties are 
known. Comparison with a known solution will show whether or not the method 
works as well as provide a benchmark on its convergence properties. Specifically, we 
consider the evolution of a constant scalar field in an otherwise empty, flat Friedmann- 
Lemaitre-Robertson- Walker (FLRW) space-time. Although extremely simple, this 
example is pedagogically useful as it serves well to elucidate the employment of the 
scheme without unnecessary complications which arise from more complex systems. As 
well, the constant field evolution is identical for both the quintessence and tachyonic 
scenario. We show how to extract various quantities for the Cauchy evolution and 
compare with the known analytic result. 

Since the source consists of only the constant scalar field, we can immediately 
set the following initial data: 



X(x) = Xo = const, ^(x) = //(x) = r7(x) = if;„(x) = 0. 



(32) 



With this prescription, the equation pairs (j27al - I27b|) and (j31al - I31b|) both take the 
form: 



1 



7 \^mn;a '^na;m\ 0, 

-^R-l [7™"V^«m]' + ^7"'^7"Vn.n^afe + 8nV{xo) = 0. 



(33a) 
(33b) 



It should be noted that in general, the consistency equations, although underdeter- 
mined, are difficult to solve even in vacuum and one must appeal to numerical tech- 
niques for solutions. 
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As well, we consider the spatially flat FLRW metric: 

10 " 
r2 
sin^ 6 

with a{to) the value of the scale factor on the initial data hyper-surface. We consider 
an inflationary scenario and prescribe: 

a\t,) =^'\ (34) 
[a2(t)],o|*=*o =2e^*°. 

At this point, all the allowable initial data has been prescribed. The Einstein-scalar 
field equations ()29all29b]) dictate that 

=Xo, 

V^(Xo) =|:, (35) 

y'(</>) ^0. 

At this point it is useful to reiterate that the above prescription is only valid ij it 
satisfies the constraint equations ()30al - l30b|) . In this example, the relevant values for 
the constraint equations are: 

^mn = 2e^"'5mn (36) 

It is a simple matter to check that these quantities indeed satisfy the constraint equa- 
tions ()33a|) and ()33b|) . 

We show both the Cauchy evolved and the analytic (e^*) square of the scale 
factor in figure |21 As expected, there is excellent agreement at small values of At. 
Deviations at larger At may be minimized by retaining more terms in the Taylor series 
()28a|) as can be seen from the vatious dotted lines in the figure (see figure caption). It 
may readily be verified that the other parameters to be evolved will display the proper 
evolution by a simple inspection of ()28al - OHej) . That is, in this simple example, all 
other parameters will retain their initial values as prescribed by (j32j) . 

4 Concluding remarks 

Both the quintessence and tachyonic scalar field, supplemented with a perfect fluid, 
were considered in the context of general relativity. The general mathematical prop- 
erties of the system, including the eigenvalue structure were studied. It is seen that 



[7afe(x)] = a^(to) 
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Scale factor squared 




1 1.2 1.4 1.6 1.8 2 

t 



o o o o o o quadratic 
...... cubic 

o o o o o o quartic 
analytic 

Figure 2: Comparison of the square of the cosmological scale factor. The dotted lines 
represent the numerically evolved Cauchy data utilising the scheme outlined in this paper 
to various orders in t — to (quadratic, cubic, quartic). The solid line represents the analytic 
result (e^*). 



the system may behave either as a two component perfect fluid or an anisotopic fluid, 
the anisotropy being due to the properties of the scalar field. Finally, in the geodesic 
coordinates, the Cauchy problem as well as the initial constraint equations have been 
derived. Assuming analyticity, the Cauchy scheme presented here is convergent. The 
scheme is iterative and may easily be executed by computer. 
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